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We investigate a nonlinear localization microscopy method based on Rabi oscillations of single 
emitters. We demonstrate the fundamental working principle of this new technique using a cryogenic 
far- field experiment in which subwavelength features smaller than A/ 10 are obtained. Using Monte 
Carlo simulations, we show the superior localization accuracy of this method under realistic condi- 
tions and a potential for higher acquisition speed or a lower number of required photons as compared 
to conventional linear schemes. The method can be adapted to other emitters than molecules and 
allows for the localization of several emitters at different distances to a single measurement pixel. 

PACS numbers: 42.50.Gy, 42.50.Nn, 42.62.Fi 



Optical microscopy has experienced many revolution- 
ary developments in the past two decades such as scan- 
ning near-field optical microscopy (SNOM) [1, 2], two- 
photon confocal microscopy, coherent antistokes Ra- 
man scattering (CARS) [3], stimulated emission deple- 
tion (STED) [4], and single-molecule localization mi- 
croscopy [5, 6]. One of the highest three-dimensional spa- 
tial resolution achieved so far has reached about 2 nm [7] 
and was based on the latter concept, where the loca- 
tions of individual molecules are determined by finding 
the centers of their diffraction-limited point-spread func- 
tions [8, 9]. The requirement for this technique is the 
distinguishability of neighboring molecules so that each 
point-spread function can be examined separately. The 
first demonstration of this concept was at low temper- 
ature [7], where the inhomogeneous distribution of fre- 
quencies in the sample provides a convenient "spectral 
identity" for each molecule. More recent approaches have 
used stochastic photo-activation or switching schemes 
for addressing single molecules [10, 11]. These room- 
temperature experiments have caused a great deal of in- 
terest for their applicability to biological systems, but 
their accuracy is limited by the number of recorded pho- 
tocycles before a single molecule undergoes photobleach- 
ing. 



In this article, we discuss a new scheme of nonlinear 
localization microscopy based on the observation of co- 
herent Rabi oscillations [12]. In addition to a first exper- 
imental demonstration, we present Monte-Carlo simula- 
tions and examine the attainable localization accuracy 
as a function of the number of recorded photons, pixel 
size, and number of Rabi oscillations. In particular, we 
compare the performance of coherent Rabi imaging mi- 
croscopy (CORIM) with standard localization methods, 
where the spatial distribution of the point-spread func- 
tion depends linearly on the signal intensity. 



Experimental localization of a single molecule via 
position-dependent Rabi oscillations 



In this work we studied the dye molecule dibenzantan- 
threne (DBATT) embedded in a n-tetradecane Shpol'skri 
matrix. The excited state in DBATT has a fluorescence 
lifetime of Ti = 9.5 ns, corresponding to a radiatively 
broadened linewidth of A/27T = 17 MHz for the zero 
phonon line (ZPL) transition. We used a CW single- 
mode dye laser (Coherent 899 Autoscan, Au < 1 MHz) to 
identify single molecules via fluorescence excitation spec- 
troscopy [13]. For all further experiments we chopped the 
beam by means of two cascaded acousto optical modula- 
tors (AOMs), achieving pulse width down to 2.9 ns and 
a signal to background ratio of more than 60 dB. Fur- 
ther details of the pulse generation scheme are explained 
in a previous publication [12]. As sketched in Fig. lb, 
a confocal microscope based on a solid-immersion lens 
enabled us to efficiently excite a molecule and detect it 
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FIG. 1. (a) Experimental setup with two cascaded acousto- 
optical modulators to obtain short laser pulses with a high sig- 
nal to background ratio, (b) The optical setup inside the cryo- 
stat with a solid-immersion lens (SIL). (c) The level scheme 
of a single molecule. 
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FIG. 2. Intensity dependence of the response from a single 
molecule in temporally integrated detection. Pulse duration 
was 4 ns, repetition rate was 700 kHz, and the intensity was 
integrated over 100 ms. Figure adapted from [12]. 

both via its Stokes-shifted fluorescence and via extinc- 
tion spectroscopy [14]. The red-shifted fluorescence was 
filtered from the excitation laser light by an optical long 
pass filter and was detected using an avalanche photodi- 
ode (APD). 

In coherent state preparation, the system can undergo 
several Rabi cycles, depending on the pulse area (A) 
which is proportional to the square root of the intensity 
and the duration of the excitation pulse. Figure 2 shows 
the experimental results for the integrated Stokes-shifted 
optical response of a single molecule excited by a series of 
4 ns long narrow-band pulses. The signal was averaged 
for 100 ms and is plotted against the excitation intensity, 
depicted on a quadratic scale. We note in passing that 
as illustrated by the dotted line in Fig. 2, a purely linear 
response as used in conventional localization techniques 
would appear as a quadratic function in this plot. 

If the focal spot of the excitation beam is scanned over 
the sample, the single emitter experiences different values 
of A at various lateral coordinates, and the emitted flu- 
orescence signal appears as several concentric rings. The 
width of each ring depends on the gradient of the inten- 
sity at that point and is the smallest at the highest slope 
of the optical point-spread function. The center of the ac- 
quired image depends on the maximum excitation inten- 
sity of the molecule and corresponds to several excitation 
and de-excitation cycles of the emitter. Figures 3a and b 
show the experimental results of such a measurement on 
a single molecule with continuous- wave (CW) and pulsed 
excitations, respectively. The occasional dark regions in 
the CW image reveal typical triplet-state blinking, while 
the full width at half-maximum (FWHM) of 370 nm in- 
dicates an almost diffraction-limited focusing spot. 

The raster-scanned image with pulsed excitation dis- 
plays higher spatial frequencies and spatial features down 
to « 40 nm. The dark stripe in the middle of the image 
was caused by a temporary spectral jump of the molecu- 



FIG. 3. Experimental raster-scanned images of a single 
molecule without (a) and with (b) pulsed excitation (4 ns 
pulses, 700 kHz repetition rate). The effectively larger point- 
spread function of the illumination spot due to the field de- 
pendence is visible. The rings around the central spot have 
widths of ~ 40 nm. 

lar resonance and provides a convenient measure for the 
background fluorescence. We also note the asymmetry of 
the images in Fig. 3, which might be caused by a cush- 
ion distortion slightly away from the optical axis of the 
solid-immersion lens [15]. 

Analytical considerations 

In Fig^ 4a, we first consider the simulated linear re- 
sponse P re d oc Iiaser of a single molecule, where P re d 
denotes the red-shifted fluorescence signal. One can lo- 
calize the molecule as is common in conventional local- 
ization schemes by applying a numerical regression to 
an assumed Gaussian intensity distribution [8, 9]. The 
highest spatial frequency in this imaging configuration is 
given by the diffraction-limited profile of the laser beam 
focus. 

We now consider coherent excitation with light pulses 
short compared to the lifetime, t p <C T\. Under this 
condition, the molecular polarization and the population 
of the ground and excited states oscillate at the Rabi 
frequency 

f2 = dE(x ,y )/h (1) 

where E(xo,yo) is the electric laser field at the position 
{#o?2/o} °f the molecule and d is the transition dipole 
moment in the direction of the field. The state at the 
end of the pulse is determined by the pulse area which 
for a rectangular pulse shape reads 

A = Qt p OC ^Iiaser- (2) 

As a result, the excited-state population at the end of the 
pulse follows a sinusoidal shape with increasing A. Ac- 
cordingly, the emission probability of a photon per pulse 
is given by 

P = (sin A/2) 2 . (3) 
Assuming a Gaussian focal shape, the electric field reads 
E( X(h yo) = E e-^ x - Xo)2 ^ y - yo)2 y^ (4) 
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where {x, y} is the position of the focal spot which is 
scanned over the x-y plane. The parameter wo gives the 
characteristic width of the Gaussian shape. Thus, finds 
the photon-emission probability to be 

P = sin 2 (/ v ^oe-[^-^) 2+ ^-^) 2 ]/^o)) ( 5 ) 

where Iq is the light intensity at the focal spot and / is 
a scaling parameter. In what follows, we use this expres- 
sion for the Monte-Carlo simulations. 

Experimental measurements such as those in Fig. 3 
might deviate from the ideal response. Usually the ob- 
served Rabi oscillations might be damped for three rea- 
sons: a) a finite pulse length and the dephasing time T2 
reduce the modulation depth below unity, b) due to the 
finite pulse length, the molecule decays and is re-excited 
within the pulse duration several times, and c) the pulse 
area fluctuations wash out the visibility of the integrated 
signal. These effects are collectively included in the fol- 
lowing equations by damping the oscillations with a field- 
dependent term a. Furthermore, we add a linear con- 
tribution b to account for the background, which scales 
roughly linearly with the optical field as experimentally 
observed in [12]. Also other molecules might spectrally 
overlap and lead to a background contribution which is 
then proportional to the excitation intensity, modeled 
by c. Moreover, the intensity-independent pixel noise 
(i.e. dark counts from a photodetector) can be denoted 
by r. With the above-mentioned parameters it is possible 
to describe the full optical response in the following form 



sin 



77- 



by/l + cl 



(6) 



where 77 covers the collection efficiency, pulse repetition 
rate, and the reduction by filters used to discriminate 
between laser and fluorescence light. Figures 4b and d 
show the contributions of the first and second terms on 
the right-hand side of Eqn. (6). 



Monte Carlo analysis of the localization accuracy 

In order to assess the promise of CORIM for local- 
ization microscopy, we performed Monte Carlo simula- 
tions using Mathematica 7.01 (Wolfram Research). For 
realistic experimental parameters, we used a statistical 
evaluation by generating images followed by a numerical 
fitting procedure. This method allows us to judge the 
localization accuracy for a variety of parameters such as 
the number of detected photons. For each data point and 
set of parameters, a series of 500 images was generated, 
by filling a grid of the defined pixel size according to the 
distribution function (see Fig. 5a) given by Eqn. (5) until 
the detected number of photons was reached. Figure 5b 
displays an example of the Monte-Carlo image. The cen- 
ter of the emitter was randomly chosen to be within the 
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FIG. 4. Simulated raster-scanned image of a single molecule 
and line scans, (a) The usual lineear optical response of a sin- 
gle molecule, (b) If the excitation light is pulsed and the pulse 
width is shorter than the singlet T\ time of the molecule, sev- 
eral rings appear, representing Rabi oscillations, c, d) Cross 
sections of (a) and (b). 
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FIG. 5. Basic derivation of pixilated images for Monte Carlo 
estimations: (a) Calculated point-spread function in the ab- 
sence of amplitude damping and detector noise (analytical), 
(b) Numerically estimated pixel distribution for 1000 detected 
photons, (c) Line cuts for pixel distribution and fitted point- 
spread function. 



inner 2x2 pixel. The pixel grid was defined such that 
the range around the molecule was larger than 1 /im by 
one pixel. The resulting image was fitted by the point- 
spread function using the Mathematica internal func- 
tion for nonlinear fitting without any further constraints 
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(NonlinearModelFit). The cross section in Fig. 5c shows 
an example. Here, experimental background noise has 
not been included so that only pixilation and shot noise 
contribute to the observed fluctuations. 

The start parameters of the simulations that fol- 
low below included the actual center of the spot, 
the height of the highest pixel as an amplitude, 
and wo of the optical point-spread function (wq = 
\ FWHM/V2 In 2 « 120 nm). The 500 Monte Carlo runs 
were unweighted averaged. The confidence range (±1 a) 
of one coordinate is displayed. Unless differently spec- 
ified, the parameters were as follows: Gaussian focus, 
FWHM=300 nm, detection of 200 photons, pixel size of 
50 nm, and a Rabi parameter of 100, leading to approxi- 
mately 3 Rabi flops (6.2 tt). These parameters are labeled 
by the dashed lines in the presented figures. 

The localization accuracy with N detected photons 
scales as 1/ a/TV, depending on the shot noise of the mea- 
surement. This behavior is observed both for linear imag- 
ing and CORIM. At the top of Fig. 6, we show CORIM 
and Gaussian simulations for different N. The bottom 
part of Fig. 6 summarizes the results and compares the 
localization error between the two cases. We find that 
the accuracy in CORIM is superior by a factor of two 
for the given parameters. Thus, a single emitter can be 
localized with higher accuracy in a shorter time, which 
might offer a crucial advantage against systematic errors 
caused by mechanical drift and jitter [16]. 
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FIG. 6. Localization accuracy for different numbers of de- 
tected photons. Gaussian fitting (upper red curve) and co- 
herent Rabi imaging microscopy (lower blue curve). Both 
curves scale as 1/y/N, as expected for pure shot-noise lim- 
ited localization performance (background noise is neglected 
in these calculations). For details of this simulation see text. 

Linear imaging requires operation below the saturation 



regime because otherwise the image of a single molecule is 
distorted with respect to the expected point-spread func- 
tion. In CORIM, on the other hand, higher intensities 
simply result in more photo cycles. In fact, the number 
of Rabi oscillations has a direct impact on the localiza- 
tion accuracy. In Fig. 7 the excitation intensity and thus 
the number of oscillations are varied. This is realized by 
changing the pulse area in Eqn. 3. For zero amplitude the 
CORIM response is equivalent to the linear imaging such 
that both curves originate at one point. As the number 
of Rabi oscillations increases, the localization accuracy 
becomes worse than for Gaussian fitting because of a vir- 
tual broadening of the recorded point-spread function. In 
other words, the same number of photons is distributed 
unspecifically on a fiat-top distribution. However, when 
the field strength is increased further, the resulting spot 
is "opened" at the center and the localization accuracy 
becomes superior to that of a Gaussian fitting. Here, 
the larger number of observed Rabi oscillations helps to 
achieve higher spatial frequencies, allowing for better fit- 
ting. The slight undulations of the localization accuracy 
versus the field strength (E(xo, yo)) shows that each time 
that the optical response opens (this occurs at odd mul- 
tipliers of 7r), the accuracy is improved. The slope of the 
curve in Fig. 7 yields a 1/ ^/E(xo, yo) behavior for the 
relative excitation intensity. 

As shown by Thompson and coworkers [9], the local- 
ization accuracy depends also on the pixel size. In Fig. 8 
we present a comparison between linear imaging and 
CORIM regarding the influence of the pixilation noise 
on the localization accuracy. We see that in both cases 
the localization procedure becomes meaningless for large 
pixels, however CORIM turns out to be more tolerant in 
this respect. 

The decisive parameters in finding the center of an 
emission spot are the total number of detected photons, 
the pixel size, and pixilation noise causing higher fre- 
quency components in the image [17]. Besides these fac- 
tors, any additional source of noise introduces extra error 
in the localization accuracy. One important omnipresent 
example stems from the background fluorescence of the 
sample and other molecules in it. In the case of pulsed 
excitation two other issues have to be kept in mind. If 
the pulse has a finite length, there is a certain probability 
that the single emitter decays and is re-excited within the 
pulse duration. Furthermore, pulse fluctuations in time 
and intensity wash out the modulation depth. 

While we have chosen to analyze Monte Carlo sim- 
ulations in this work, a variety of other measures can 
be used for evaluating the localization accuracy in mi- 
croscopy [9]. For example, one could start with the ex- 
pected point-spread function (see Fig. 5a), add different 
kinds of noise at each pixel and analyze the localization 
error. A similar strategy corresponds to the characteri- 
zation of the Cramer-Rao lower bound by calculating the 
Fisher information contents [18]. However, regardless of 
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FIG. 7. Localization accuracy for Gaussian localization (red 
curve) and for coherent Rabi imaging microscopy (blue curve). 
With increasing pulse strength, the Gaussian point-spread 
function widens first and reduces the localization accuracy 
by increasing the effective spot size. At an intensity of about 
27T, the localization gets better than linear imaging. The oscil- 
lations for higher field strength are explained by opening and 
closing the inner spot of the point spread function. See text 
for details on the excitation strength units and the details of 
the simulation. 
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FIG. 8. Localization accuracy for different pixel sizes. The 
coherent Rabi imaging microscopy (blue curve) has a higher 
stability against fitting. Note that background noise is ne- 
glected such that only pixilation noise is present. For details 
of this simulation see text. 
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FIG. 9. (a) Simulated image of two molecules located at a 
distance of 20 nm experiencing the same effective pulse area 
with the same illumination strength. The concentric rings 
with sub- wavelength features can be easily differentiated, (b) 
A line cut through the central of the point-spread function. 



the exact procedure, it is apparent that higher spatial 
frequencies in the overall point-spread function provide 
more information and improve the localization accuracy. 
In particular, the largest slope of the point-spread func- 
tion contributes most to the localization accuracy [8]. 



Resolving several emitters 

Having discussed the accuracy in the localization of a 
single molecule, we now explore the potential of CORIM 
for resolving several close-lying emitters. Since in the ex- 
periment we only collect red-shifted incoherent photons, 
we simply add the amplitudes and do not need to include 
interferences [19]. 

Figure 9 displays an example of two emitters sepa- 
rated by 20 nm along the x-axis. On the separation axis 
it is apparent that one maximum of the response coin- 
cides with the minimum of the other emitter. The situa- 
tion is analogous to magnetic resonance imaging (MRI) 
or super-resolution microscopy in inhomogeneous electric 
fields [7], where externally imposed static field gradients 
lead to position-dependent spectral shifts. Here the field 
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gradient is provided by an inhomogeneous light intensity 
distribution. If we consider several concentric rings, we 
can assign a conservative bound to the localization accu- 
racy by the resulting width of the rings. 

Let us start by considering emitters with negligible 
transition frequency differences. As a concrete exam- 
ple, we assume four close- lying emitters as depicted in 
Fig. 10a and recorde the Rabi response at a central pixel 
labeled by a green cross. At this point the emitters feel 
different effective pulse areas and thus undergo different 
Rabi frequencies, leading to a complex beating pattern 
that results if one plots the fluorescence signal as a func- 
tion of field strength (see Fig. 10b). However, Fig. 10c 
reveals that the Fourier transform of this signal clearly 
identifies four Rabi frequencies corresponding to the four 
molecules separated by a fraction of the wavelength. 

If the transition frequencies of the emitters are 
slightly different (as e.g. in the inhomogeneous band of 
molecules), it is still possible to address them simulta- 
neously by using an optical pulse that is short enough 
to cover their frequency differences. If the emitters show 
a considerable inhomogeneous broadening, the frequency 
detuning A contributes to the effective Rabi frequency 
according to i? e ff = \J Q 2 + A 2 . If we now scan the fo- 
cus across the sample, the Rabi response differs not only 
because of the effective excitation strength, but also be- 
cause of the different spectral response. As a result, the 
asymmetry of the point-spread function is larger than in 
the earlier case and lateral shifts in the order of 1 nm can 
be resolved. This also represents a convenient method to 
locate coupled molecules as it has been observed in [7]. 



Outlook 



We have introduced a method to localize emitters 
based on coherent Rabi oscillations. The technique in- 
troduces higher spatial frequencies in the recorded im- 
ages, allowing for a higher localization accuracy with 
fewer photons. It is possible to extend these experiments 
to ambient conditions if the coherent features of single 
molecules are explored further [20]. Moreover, a single 
emitter can be used to sample and characterize an opti- 
cal pulse at nanoscopic scales. 

The localization technique discussed here can be read- 
ily generalized to other coherently driven processes and 
emitters. For example, one can apply it to ions in a 
trap where numerous Rabi oscillations have been demon- 
strated with a very high fidelity [21]. Furthermore, 
CORIM can be applied to emitters such as NV- centers in 
diamond with microwave transitions of the ground state. 
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FIG. 10. (a) Simulated image of four molecules (red dots), 
located at a small distance from the measurement point. The 
field dependence measured on a central pixel (green cross) is 
displayed in (b) . c) The fourier transform of the Rabi oscilla- 
tions reveal that although measured only at one single pixel, 
all four emitters can be determined. If the inhomogeneous 
broadening of the Rabi response is negligible, the distance of 
all emitters to the measurement spot can be determined. 
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